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ABSTRACT 

Context. A recent JHR^ study of several grand-design spiral galaxies, including NGC 2997, shows a bimodal distribution of their 
system of star clusters and star forming complexes in colour-magnitude and colour-colour diagrams. In a comparison with stellar 
population models including gas, the {J - H) vs (H - K^) diagram reveals that embedded clusters, still immersed in their parental 
clouds of gas and dust, generally have a redder (H - K^) colour than older clusters, whose gas and dust have already been ejected. 
This bimodal behaviour is also evident in the colour-magnitude diagram M^ vs (J - K^), where the brightest clusters split into two 
sequences separating younger from older clusters. In addition, the reddening-free index Q^ ={H- K^) - 0.884 (7 - H) has been shown 
to correlate with age for the young clusters and thus provided an effective way to differentiate the embedded clusters from the older 
ones. 

Aims. We aim to study the behaviour of these photometric indices for star cluster systems in the Local Group. In particular, we in- 
vestigate the effectiveness of the 2d index in sorting out clusters of different ages at their early evolutionary stages. In addition, the 
whole set of homogeneous measurements will serve as a template for analyses of the populations belonging to distant galaxies that 
are unresolved clusters or complexes. 

Methods. Surface photometry was carried out for 2MASS images of populous clusters younger than ~ 100 Myr whose ages were 
available. The integrated magnitude and colours were measured to a limiting radius and combined to generate the photometric dia- 
grams. Some clusters, particularly the embedded ones, were studied for the first time using this method. 

Results. The integrated magnitudes and colours extracted from the surface photometry of the most populous clusters/complexes in the 
Local Group shows the expected bimodal distribution in the colour-colour and colour-magnitude diagrams. In particular, we confirm 
the index Q^ as a powerful tool for distinguishing clusters younger than about 7 Myr from older clusters. 

Key words, galaxies: star clusters: general - infrared: galaxies - techniques: photometric 
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. 1. Introduction clusters with Qd > 0.1 concentrate around (7 - /i:)x! 1.8, with 

'^ significant scatter due to variable extinction. They form a well- 

O The study of 55 clusters that present strong Br^ emission ^^^^^^ branch that is clearly separated from the bright older 

rn : (Br^ - /: < -O.l) m the grand-design spiral galaxy NGC 2997, ^^^^^^^ ^^^^ concentrate around (J - K)^l.2. This two-branch 

T^^ , observed with HAWK-I at VLT, has found a relationship ^^^^^^.^ j^ common to ten grand-design galaxies analysed by 

;> ■ between the reddening-free index ga and the Br^ index Grosb0l & Dottorii dfoTl and even for NGC 7424 whose clus- 

■ (| Grosb0l&Dottori| | 2O12D : (Br^ - K) - -0.48 Qa - 0. 17, in- j^j. population is 3 magnitudes fainter than that of NGC 2997, as 

dicating that Qa ^ (H - K) - 0.844 {J - H) represents a good ^een in their figures 8, 9, and B.3. If the extinction varied only 

H age indicator for Qa > 0. 1 which corresponds to a^es < 7 Myr ^j^^jy ^ j^^ age, one would not expect a separation between the 

y. according to Starburst99 C Leithereij | 1999[ [ Vazquez & Leithereij y^^^g ^^^ ^j^ branches in their diagrams. Thus, the gap sug- 

120051 hereafter SB99) evolutionary tracks. Population models gg^^^ ^ ^^^^^^ change in dust properties and/or distribution that 

generated with SB99 show nebular emission for the youngest j^^^ ^^ ^ ^.^pj^j reduction of the extinction. 

clusters. The emission decreases with time abruptly around 7 

Myr, depending on model parameter s, as discussed in Section [H Studies of our Galaxy, the Large Magellanic Cloud (LMC), 

where errors to this value are also set. l Grosb0l. D ottori & Gredell and the Small Magellanic Cloud (SMC) show that the interplay 

(|2006|) previously calibrated the equivalent width of Br^, or between stellar evolution and dust properties in the intracluster 

equivalently, the Bry index (Br^ - K), as an age indicator based medium is a dynamic one: stellar clusters are bom embedded 

on the ISAAC/VLT K-band spectra and SB99 evolutionary within giant molecular clouds. During their formation and early 

tracks. The Br^ equivalent width is antic on^elated with the emit- evolution, they are often only visible at i nfrared wavelengths , 

ting region age, as previously modelled for other H emission since they are heavily obscured by dust (i Lada & Ladall2003h : 

lines in the visible (Dottori 1981; Alonso-He n'ero et al. 1996). moreover, the embedded cluster birthrate ex ceeds that of the vis- 

lGrosb0l & Dottorii (1201 2b also found that in the (7 - K) vs Mk ible open cluste rs by an order of magnitude (iLada & L ada 20031). 

colour-magnitude diagram (CMD) of NGC2997, the youngest The findings of lGrosb0l & Dottorii (1201 2l) could help in under- 
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standing the behaviour of the intracluster medium, which could 
be associated with the first supemovae explosions and stellar 
winds from O-type stars, when dust was expell ed from the clus- 
ter environment JB astian & Goodwin 2006; Go odwin & BastianI 
I2006t iLada & Lada. ,2003; .Whitworth 1979). The increase in 
{H — A's) for younger clusters is also a consequence of the pres- 
ence of hot stars because they ionise hydrogen, producing con- 
tinuum emission and the Br^ line, and also heat dust, both con- 
tributing to the flux in K. 

We checked the capacity of the near-infrared (NIR) diagrams 
to differentiate clusters younger than 7 Myr from the older ones 
in a sample of clusters from the Milky Way, the Magellanic 
Clouds (MCs) and two clusters from M 3 1 and M 33, all of which 
have well-determined ages. In the next section, we present the 
NIR photometry of the Local Group clusters obtained from the 
2MASS images. In Section 3, we compare our mul ti-diaphragm 
phot ometry of LMC-hSMC clusters with that of iPessev et al.l 
(l2006i hereafter P06). The age limit estimate and the NIR colour 
sensitivity on parameters/ingredients of stellar populations evo- 
lution models are presented in Section 4. In Section 5, we dis- 
cuss the observational diagnostic diagrams. In Section 6, we 
show how the gd index can be used as an age diagnostic tool. 
In Section 7, we discuss the results and in Section 8, we present 
our conclusions. 







log[age(yr)] 

Fig. 1. Absolute K magnitude versus age for the Local Group 
clusters in our sample. 



K magnitude with age for the Local Group clusters, characteris- 
ing the sample in terms of intrinsec properties. 



2. Near-infrared surface photometry of populous 
young star clusters 

The sample listed in Table [T] is based on the selection by 
iPortegies Zwart et al.l (l2010l) . complemented with the Milky 
Way embedded clusters RCW 38, NGC 3576 and M 1 7, which 
have masses higher than 500 Mq dLada & Ladall2003h . and the 
LMC clusters NGC 1805, 1984, 201 1, 2156, 2159 and 2 172, 
which were studied by P06. .Portegies Zwart et al.l (1201 Ol) re- 
stricted their sample to well-studied populous clusters with 
known ages; in addition to the MCs and Milky Way clusters, 
NGC 604 in M33 and vdbO in M31 were also analysed. All 
clusters have ages up to 100 Myr. Such a sample was considered 
suitable to derive the integrated light properties using 2MASS 
mosaic images built from relatively shallow exposures. The co- 
ordinates, age, mass, reddening and distance adopted for each 
cluster are given in Table [T] 

The JHK^ images of the 42 star clusters were downloaded 
from the 2MASS database using the software MontagqJ, which 
assembles single images (512 x 1024pix; l"/pix) into custom 
mosaics preserving the spatial and calibration fidelity of the in- 
put images. A zeropoint keyword determined by the average of 
zeropoints provided in all single images was appended to each 
mosaic file for each band. Typically, the mosaics cover sky areas 
from 2° X 2° to 12' x 12' depending on the cluster size. 

The clusters' centres were determined visually in the images, 
and the surface brightness was evaluated in the annular regions 
around this centre. Fixed width annuli were employed for each 
cluster, but they varied from 2" to 10" depending on the cluster 
size. The sky background (Ft,g + o-[,g, in DN/pix) was determined 
for each band using the whole mosaiced frame by means of an 
algorithm that evaluates the mode of the sky distribution and re- 
jects pixels contaminated by stars. 

The cluster integrated colours and their uncertainties, the 
adopted aperture radius and the calculated absolute K magnitude 
are given in Table[T] Fig.[T]shows the distribution of the absolute 



3. Comparison: this work vs Pessev et al. 2006 

Fig.|2]shows the measured JHKs integrated magnitudes of MCs' 
star clusters in common with those in the P06 sample. In the 
present study, we defined the area covered by a cluster as limited 
by the radius for which the magnitude difference from consec- 
utive diaphragms is approximately the same as the magnitude 
error. The photometry performed by P06, although also based 
on 2MASS images, followed a different approach: both the PSF 
photometry (for the resolved component) and the surface pho- 
tometry (for the unresolved component) were employed, and 
probable non-members were excluded from the final set of aper- 
ture measurements. The procedure caiTied out by P06 to exclude 
probable non-member bright stars from the integrated photome- 
try was based on a comparison between the individual stars mag- 
nitudes derived from PSF photometry and their expected locii 
in a CMD considering the cluster age. Discrepancies between 
them indicated a probable non-member, which was then sub- 
tracted from the cluster integrated light. In our measurements, 
we did not exclude any star because our intention was to com- 
pare the clusters' integrated light with those obtained for clus- 
ters/complexes observed in distant galaxies. While P06 included 
different error sources, even the eiTor introduced by background 
stochastic fluctuations, in the present analysis, only the photo- 
metric errors were considered. To quantify these eiTors, we used 
for each band the total noise variance as given in the Explanatory 
Supplement to the 2MASS All Sky Data Releasfl i.e., cri,a„^ = 

where n is the number 
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n X 0-2^ X [(4x1.72) -Hi. 8x10- 



'«], 



of pixels in the aperture and (Tbg is the measured background 
noise (in DN/pix). The magnitude error in a given aperture was 
then obtained by o-,„„g = 1.0857 x o-ba„d/Fha„d, where Ftand is 
the background subtracted flux in DN. The calibrated integrated 
magnitude in each band was magcai - zeropoint - 2.5 log Fhand- 
The colour errors were determined by propagating the integrated 
magnitude errors from the relevant bands. Because P06 yields 
integrated magnitudes for several diaphragms, we selected the 



' Montage v3.3 is available in the NASA/IPAC Infrared Science 
Archive ( http://irsa.ipac.caltech.edU/applications/2MASS/j 



http://www.ipac.caltech.edu/2mass/releases/allsky/doc/sec6_8a.html 
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Fig. 2. Comparison of our integrated near-IR magnitudes of MC 
clusters with those from P06. J, H and Ks magnitudes are plotted 
as different symbols and colours. The straight continuous line is 
a one-to-one relation. Line segments connect the three magni- 
tudes for exemplary clusters. 
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Fig. 4. Colour-colour diagram of the outlier clusters indicated in 
Fig-E] P06 measurements are indicated by filled green symbols. 
Large open red symbols represent the effect of adopting differ- 
ent cluster centres. Small open blue symbols show the effect of 
different sky backgrounds. 
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Fig. 3. Comparison of our integrated near-IR colours of MC 
clusters with those from P06. (7 - H) and (7 - K^) colours 
are plotted as blue diamonds and red squares, respectively. The 
straight continuous line is a one-to-one relationship. Line seg- 
ments connect colours of cluster outliers. 



one that matches our adopted limiting radius for the following 
analysis. 

The integrated magnitudes comparison is shown in Fig. |2] 
There is a general agreement between both studies for all three 
bands. The continuous line is a one-to-one relation. Some clus- 
ters, taken as examples to illustrate the colours behaviour, are 
indicated by line segments connecting the three bands. Line 
segments parallel to the one-to-one relation reveal clusters that 
have identical colours in P06 and our study. Crooked or slanted 
line segments yield different colours in P06 and our study. The 
colours comparison is shown in Fig.|3]and indicates clusters with 
dissimilar colours compared to those from P06. 

The effect of the different adopted centres and background 
corrections were investigated for these clusters to search for an 
explanation for the discrepancies found. Fig. |4]pi"esents the re- 
sults of the changing centre and the background correction for 
the aforementioned clusters on the (7 - H) vs (H - K^) diagram. 



The adopted cluster centre plus four positions shifted from this 
centre by 10" towards the E, W, N and S, respectively, were 
considered, and the colours were determined for each of centre. 
The sky estimates were evaluated in five different fields, defined 
as square areas of a typical size 2 times larger than the cluster 
adopted limit (see Tab. [TJ, and the colours determined for the 
clusters maintained their original adopted centres. The colours 
derived from the original sky estimate using the whole mosaiced 
frame is also plotted, as well as the colours from P06 for these 
four clusters, where we plotted their colours for the diaphragm 
that matches our adopted cluster limit. The effect of the moder- 
ate centre shifts and the different sky corrections appear not to 
account for the discrepancy of the cluster colours in comparison 
with the P06 measurements. Particularly, NGC2172 presents a 
huge difference in both colours. NGC2156, on the other hand, 
in spite the larger photometric errors, present colours compati- 
ble with P06. NGC 1818 and NGC 1847 have a disagreement in 
{H - K,). 

Fig. |5] shows the (J - H) vs (H - Ks) diagram for the same 
clusters but with integrated colours calculated from the 2MASS 
point source catalogue (PSC); i.e., the individual star flux for all 
stars within the adopted cluster limit was added up. The colours 
resulting from the different centring (10" shifts) were also plot- 
ted. The different centres do not affect the cluster colours sig- 
nificantly. Because P06 did not provide a list of members and 
the star distribution in the NIR CMDs of these clusters does not 
define well the expected evolutionary sequences, we did not at- 
tempt to subtract the integrated flux from the flux of the possible 
non-members, which prevented us of making a direct compari- 
son with the P06 photometry. 

It is worth noticing, however, that the four clusters men- 
tioned with colour discrepancies have no apparent sign of 
gas. We searched for stars that c ould be respons i ble fo r the 
colour differences. In accord to IWorthev & Leel (l201ll) and 
IStraizvs & Lazauskaitd (l2009h there are no stars with colours 
blue enough to produce a.{J - H) colour like that of NGC 2172 
quoted by P06. Our integrated colours calculated from surface 
photometry (Fig.|4]l and individual star fluxes (Fig.|5) agree, but 
diverge from that of P06. As NGC 2172 colours would fall in 
a rather peculiar position in the colour-colour diagram and not 
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Fig. 5. Colour-colour diagram of the outlier clusters indicated in 
Fig.|4]based on the integrated flux of individual stars within the 
adopted cluster limit. The symbols are the same as in Fig.|4] The 
influence of the diff'erent centres is shown for each cluster. 



being fit by any model (see Fig.|9]l, this indicates that P06 pho- 
tometry may be incorrect for NGC 2172. Because the differences 
between our photometry of NGC 1818 and NGC 1847 and that 
of P06 are not as large as for NGC 2172, their colour discrep- 
ancies may be explained by the exclusion of bright non-member 
star(s) in P06 photometry. NGC 2156 has a smaller colour dis- 
crepancy which also can be attributed to non-member star(s) be- 
ing removed in P06 photometry but not in ours photometry. 



Fig. 6. Equivalent width of Biy as a function of age according 
to SB99 models. The base model parameters/input are defined 
by continuous star formation (CSP), Kroupa IMF with upper 
mass m„=lOOM0, metallicity Z=0.02 and Padova tracks (con- 
tinuous black line). In the following models, all parameters/input 
of the base model are kept, except for one variable replace- 
ment: simple stellar populations (SSP) (dotted red line); Geneva 
tracks (dot-dashed blue line); upper mass m„=120MQ (dashed 
cyan line); upper mass m„=3OM0 (dashed red line); metallic- 
ity Z=0.008 (dot-dashed green Une); metaflicity Z=0.05 (dot- 
dashed magenta line); Kroupa IMF with upper slope x=2.0 (dot- 
ted blue line); Kroupa IMF with upper slope x=2.6 (dotted green 
fine). 



4. Evolution of stellar populations NIR colours: 
model predictions 

On theorethical grounds, it can be demonstrated that the drop 
of Bry emission, associated with the evolution of hot stars that 
are not capable anymore of ionising the gas in the intraclus- 
ter medium, occurs at an age of ~7 Myr. With this aim. Fig. |6] 
presents the behaviour of the Br^ equivalent width with the age 
according to SB99 models (v6.0.4). We used as base model the 
one with con tinuous star formation population (CSP), standard 
iKroupa (2001) IMF (x=2.3 for masses above 0.5 Mq and x=1.3 
for masses below 0.5 Mq) with upper stell ar mass m„=100Mo , 
solar metallicity (Z=0.02) and Padova jBerteUi et al.1 11994 
iGirardi et al.ll2000i) stellar evolutionary tracks. Other models for 
which one parameter or ingredient is changed in the base model 
are represented in Fig. |6] The whole set of models entails the 
effect of several variables on the age at which the Bry drops, ac- 
cording to SB99 models. We explored a range of properties that 
represent a realistic approach to our problem. 

Based on the models considered in Fig.|6] we computed, for 
each model, the age at which the Bry equivalent width falls by 
half of its maximum value. Then the average age and its standard 
deviation were estimated as 7.5 + 0.9 Myr and considered as the 
age limit separating younger populations with nebular emission 
from older populations without nebular emission. 

SB99 models also allowed us to investigate the behaviour of 
the reddening-free index 2d with age and its sensitivity to the 
aforementioned parameters (Fig.|7). The age limit and its uncer- 
tainty as calculated above was used as a constraint for the range 
of 2d representing the transition between young and old clus- 
ters. The higher density of models intersecting the age transition 
region lies around Q^ = 0.1, consistent with iGrosb0l & Dottoril 
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Fig. 7. Reddening-free parameter as a function of age according 
to SB99 models. The model parameters/input are defined as in 
Fig. |6l The grey area represents the age (and uncertainty) esti- 
mate for the evolutionary transition between young clusters with 
nebular emission and older clusters without nebular emission. 



( 1201 2h findings. However, as a cautionary note, it can be seen 
that in this region 2d as low as -0.15 or as high as -1-0.15 are ac- 
cessible to specific models. Lower values of Q^ favour the base 
model with the SSP instead of CSP as the star formation mode 
and a low upper stellar mass. Higher values of gd favour the base 
model (CSP, m„ = 100 Mq). 

Aiming at a comparison with the observations, the parameter 
space covered by different models was further investig ated in the 
(J - H) vs {H - K,) diagram (Fig. ID. SSP models bv lMarastonI 
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Fig. 8. Colour evolution of stellar populations up to 100 Myr according to different models and ingredients. Fixed parameters are 
indicated at the top while variable ones at the bottom, (a) Comparison between SSPs built with Padova tracks and SSPs from 
iMarastonI (l2005l) for different metallicities. (b) SSPs with Padova tracks including nebular emission compared to those without 
nebular emission, (c) Same as in (b) for CSPs. (d) CSP models built with a Kroupa IMF and variable upper stellar mass limit, 
(e) CSP models built with different metallicities. (f) CSP models built with different slopes for the IMF above 0.5 Mq. Age (Myr) 
is indicated at selected evolutionary stages in panels (b) and (c). 



(I2005L hereafter M05) are shown for metallicities Z=0.009. 0.05 would fit complexes in the Antennae (after reddening correc- 

and 0.02 (solar) and a Kroupa ( 2001 ) I MF and compared to SSPs tion) better than SB99 models, which they attribute to the differ- 

built with updated Padova tracks (Marig o et al.l 120 08) for the ent treatment of the red supergiant phase. We found that Padova 

same metalli cities (Fig. [8h ). Thes e models do not include nebu- SSPs show a trend with metallicity opposite to those of M05 

lar emission. iBastian et al.l (l2006h showed that the M05 models SSPs (see the solar and above solar paths in Fig. |8^) but none 
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of them reach extreme red colours. In any case, M05 models do 
produce redder colours for the oldest SSPs compared to Padova 
SSPs, but their SSP colours are similar for populations younger 
than 7 Myr. Fig. [Sj) presents the effect of nebular emission on 
Padova SSP colours. Our base SB99 model, i.e., CSPs built with 
Padova tracks, solar metallicity, Kroupa's IMF with stellar mass 
upper limit IOOMq and nebular emission is compared to that 
without nebular emission in Fig.|8};. Clearly the effect of nebular 
emission is to increase (H - K^) and (J - H) for younger pop- 
ulations. We also investigated the base model colour sensitivity 
to the upper mass limit (Fig.|8}l), to the metallicity (Fig. [8^), and 
to the IMF slope above O.SMq (Fig.[8f). {H - K^) increases for 
the younger populations when the upper mass limit increases, 
the metallicity decreases and the IMF slope decreases. For the 
range of parameters considered, the IMF slope seems to have 
the smaller effect on the colours, but they are strongly affected 
by the upper mass limit. In general, plausible variations of the 
parameters/ingredients investigated are not capable of matching 
all Local Group cluster colours without dust extinction (see next 
Sect.). 

5. Local Group clusters: the age effect on infrared 
colour indices 

Our sample is characterised in Figs. |9]and[T0] where the mea- 
sured NIR fluxes for the 42 clusters are complemented by the 
ages found in the literature (see TablelTJ. The age provided come 
from CMD fittings of isochrones derived from different stellar 
evolutionary models and the quoted uncertainties from the lit- 
erature (large in most cases) should account for discrepancies 
resulting from the use of the different methods of age determi- 
nation, minimizing the effect o f their heterogeneity. We dete cted 
an error in the age quoted by iP ortegies Zwart et al.l (l20I0b for 
Westerlund 1, that should be 5 Myr instead of 3.5 Myr. 

The observed colour-col our diagr am (Fig, p} corr oborates 
the explanation provided by lGrosb0l & Dottoril (l20I2b for the 
distribution of the clusters in grand-design spiral galaxies. The 
clusters can be separated into two main regions in the diagram. 
The redder group is younger than 7 Myr and associated with high 
extinction, and the bluer group is older and associated with rela- 
tively low extinction. Figure |9]clearly shows that our sample of 
clusters younger than 7 Myr fills regions of high extinction that 
are not reached by their older counterparts. 

In Fig. |9l the colours of the Padova single-burst stellar pop- 
ulations (SSP) with ages between 1 Myr and 100 Myr and the 
SB99 CSP with ages between xiOMyr and 100 Myr are shown. 
Dust and nebular emission are included in the SB99 models. 
Both models were built using as ingredients a standard Kroupa's 
IMF and solar metallicity. The different extinction models are 
indicated for a visual extinction of Ay - 5. The standard red- 
dening vector, the so-called screen m odel, accounts for the fore- 
ground extinction (Indebeto uw et al.l l2005). while the reddening 
vector characterised by the dusty model accou nts for intraclus- 
ter dust, where the dust and stars are mixed (.Israel et al.lll998t 
IWitt. Thronson & Capuanoiri992h . The observed colours of the 
reddest cluster in our sample. Quintuplet, were connected to the 
colours of the SB99 CSP, with an age of 4 Myr, the cluster age 
(see Table[TJ. Quintplets' extinction corresponds to Ay ~ 16 and 
closely follows the screen dust model, as observed by comparing 
the cluster extinction slope with that of the extinction models. A 
similar argument applied to the other two clusters close to the 
Galactic centre. Arches and [DBS2003] 179, leads to the same 
conclusion: because they are younger than 7 Myr, linking their 
observed colours to their intrinsic colours (CSP models between 
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Fig. 9. Sample observed colour-colour diagram with different 
symbols indicating age ranges: f(Myr)< 7 (blue circle), 7 < 
/(Myr)< 20 (green square), 20 < f(Myr)< 100 (red diamond). 
The three young clusters near the Galactic centre and M 42 are 
labelled. Continuous and dashed lines represent the colours of 
the Padova SSPs and the SB99 CSP, respectively. Model ages 
range from very early phases to 100 Myr, indicated for the CSP 
models. The magenta dot-dashed line connects the observed and 
the intrinsic colours of Quintuplet. Blue dotted lines show the 
extinction yielded by Ay = 5 for both a standard screen redden- 
ing vector and a dusty cluster medium. 



OMyr and 10 Myr) yields extinction vector slopes more akin to 
the screen model than to the dusty model. 

Many clusters are compatible with both the screen and dusty 
models. The SB99 CSP models were generated with a stellar 
upper mass limit of m„ = 100 M©, which, in the colour-colour 
diagram, are suitable for reproducing the colours for most of the 
sample. NGC889, NGC864 and NGC2156 are better charac- 
terised by models with m„ = 20 - 30 Mq and a low extinc- 
tion, while M42 is only characterised by a young population 
(0-2 Myr), but its colours are uncompatible with normal ex- 
tinction. 

The distribution of the clusters in the {J-K) vs Mk CMD, to- 
gether with the Qi vs Mk diagram, is shown in Fig.[TOl where the 
grey rectangle indicates values of Qd corresponding to the tran- 
sition between young and old clusters. Note that the most proba- 
ble value is Q^ =0.1 and its distribution is not symmetric around 
that value. For the (J - H) vs {H - K^) diagram, Fig.fTOlreveals 
that most of the youngest clusters (blue circles) have a tendency 
to be located in specific regions of both diagrams. In particular, 
10 out of 19 clusters with ages less than 7 Myr have a reddening- 
free index of Qi > 0.1. This proportion does not change if we 
use 7.5 Myr as the age limit. The three very young bright clus- 
ters with gd < -0.3, Arches, Quintuplet and [DBS2003] 179, 
are objects close enough to the Galactic centre to have their evo- 
lution influenced by the strong tidal field, with their gas and dust 
akeady being stripped in the early phases. They appear very red 
because of the strong foreground extinction instead of their intr- 
acluster gas/dust content. If these clusters are disregarded, then 
w 60% of the youngest clusters have a reddening-free index of 
Gd>0.1. 

Because less massive clusters colours may be affected by 
stochastic effects arising from small numbers of bright stars, a 
new colour-colour diagram was built (see Fig. \TT\ by binning 
in age the cluster colours, computed from the integrated flux 
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Fig. 10. Left: Absolute Kg band magnitudes versus the 
reddening-free index Qd for 42 clusters in the Local Group with 
ages up to lOOMyr; The grey bar indicates the range of 2d val- 
ues corresponding to the transition between young and old clus- 
ters, according to Fig. |7] Right: Colour-magnitude diagram for 
the same clusters. 
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Fig. 11. Colour-colour diagram for age binned clusters with bin 
log f=0.2. The age (log t) is indicated. Padova SSP and SB99 
CSP models are superimposed. Dust effect according to the dusty 
and the screen models are represented as in Figj9] 

in JHK and consequently these are flux-weigthed colours. This 
procedure guarantees that intrinsically fainter clusters will have 
low weight. To compute the final colours we did not weight the 
individual cluster colours by their signal to noise or any other 
parameter. The age bin is 0.2 in log t. There is a clear separation 
among the binned clusters indicating two sequences in that the 
younger populations are redder in (H - K^). 



6. Age diagnostic based on Q^ 

The age distribution of clusters up to lOOMyr can be divided 
into clusters older (O) and younger (Y) than 7 Myr according to 
Qd < 0.1 and Qd > 0.1, respectively. Applying this to our ob- 
served sample of Local Group clusters, for Qd < 0.1, we detect 
21 true old clusters (Ot) and 6 (excluding the three clusters near 
the Galactic centre) false old clusters (Of), i.e., young clusters 
with 2d < 0.1. Similarly, for Qd > 0.1, we detect 10 true young 



clusters (Ft) and 2 false young clusters (Ff), i.e., old clusters with 
Qd > 0.1. Consequently, 78% of the clusters having Qd < 0.1 
are true old clusters (Ot/(Oi + Of)), and 22% are misclassified 
clusters younger than 7 Myr. In a similar way, 83% of the clus- 
ters having Qd > 0.1 are true young clusters (Ft/(Ft + Yf)), and 
17% are misclassified clusters older than 7 Myr. The statistics 
also shows that, by using Qd = 0. 1 as a discriminator, we recover 
^ 60% (YJ(Yt + Of)) from the total sample of young clusters and 
» 90% (Ot/(Ot + Yf)) from the total sample of old clusters. 

Because we selected the most massive clusters and star form- 
ing complexes of the Local Group galaxies, the sample is, on av- 
erage, representative of the brightest (youngest) clusters in sim- 
ilar environments. Consequently, beyond the qualitative agree- 
ment of our sample NIR diagrams with those of distant galaxies 
( lGrosb0l & Dottoril2012 ). the Qd index provides a reliable quan- 
titative measure of the proportion of clusters/complexes younger 
than 7 Myr in galaxies. A large fraction (83%) of the clusters 
with Qd > 0. 1 are genuinely younger than 7 Myr and a similarly 
large fraction (78%) of the clusters with Qd < 0.1 are in the age 
range 7- 100 Myr. 

7. Discussion 

Our results validate the use of the Qd index to differentiate 
star clusters in their early stages of evolution. The sample em- 
ployed in the present study is not complete, but by selecting 
the most massive clusters/complexes with available ages in the 
Local Group, we hope to have gathered an unbiased sample. One 
difference from the cluster samples in Grosb0l & Dottori (201^ 
is the lower cluster mass limit in our sample (« 500 Mq), which 
may introduce some bias due to the integrated colour scattering 
produced by a few heavy OB stars. Another issue is that, be- 
ing lower metallicity galaxies, the MC clusters may have lower 
metallicity, in average, than those in the Milky Way and M31, 
making the sample not completely homogeneous on this aspect. 

The Qd vs Mk diagram (Fig.fTOt clearly shows a segregation 
of clusters by age and led us to give an estimate of the proportion 
of clusters younger than 7 Myr over all clusters with Qd > 0.1, 
which should be expected in galaxies outside the Local Group 
if the processes involved in cluster formation and evolution are 
similar. As Fig. [TO] shows, there is very little contamination of 
clusters older than 7 Myr among those with Qd > 0.1. It is fea- 
sible to estimate the size of the population of clusters younger 
than 7 Myr for a distant galaxy observed in JHK^. From the 
measured colours of the star clusters and the star forming re- 
gions and the calculated Qd index, one can count the number of 
clusters with Qd > 0.1, which are x; 60% of the total number 
of objects with 7 Myr. Consequently, the total number of these 
objects can be obtained. A similar analysis for Qd < 0.1 would 
lead to an estimate of the population of older clusters. Because of 
their peculiar evolutionary characteristics, the clusters near the 
Galactic centre were excluded from the above statistics, and the 
object counts in distant galaxies should exclude those clusters 
too close to their centres. Nuc lear clusters were also excluded in 
the lGrosb0l & Dottoril (1201 2l) study of cluster systems in exter- 
nal galaxies. 

We showed that the clusters near the Galactic centre in our 
sample, i.e.. Arches, Quintuplet and [DBS2003] 179, have high 
extinctions associated mainly with foreground reddening (as 
opposed to intracluster dust), which nearly follows a standard 
screen reddening vector (Fig.|9]l- They are clusters younger than 
7 Myr and their extreme NIR colours are a consequence of the 
heavy extinction towards the Galactic centre. The stellar popu- 
lation at these early ages cannot produce extremely red colours. 
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as demonstrated by the locus of the evolutionary models plotted 
inFig.ia 

In summary, if Qd > 0.1, then 60% of the clusters must 
be genuinely younger than 7 Myr, with very few older clusters 
showing up in this region (Fig. |9]l. The remaining 40% of the 
clusters, with an index of Qd < 0.1, can be explained by their 
association to extinctions consistent with a standard screen red- 
dening vector Such very young, dust depleted clusters, as in the 
cases of the three clusters near the Galactic centre in our sample, 
can occur as a consequence of extreme environmental conditions 
such as strong tidal fields. 

8. Conclusions 

We analysed the NIR colour indexes of 42 clusters in the Local 
Group with ages determined by other authors to test the capac- 
ity of the reddening-free colour index Q^ to distinguish clusters 
younger than 7 Myr from the older clusters. The index failed to 
classify 6 out of the 16 clusters younger than 7 Myr and mis- 
classified 2 out of the 23 clusters older than 7 Myr Because the 
reddening-free index verifies Q^ > 0.1 for a broad variety of 
galaxies, cluster masses and environments, we conclude that it 
is most likely associated more with the dust distribution with re- 
gard to the cluster stars than with the metallicity or the IMF. 

The sudden change of the index is most likely linked 
to the expulsion of the dust from the cluster due to stellar 
evolution, mainly SN explosions, characterising th e transition 
from the mix of stars and dust inside the cluster (llsrael et al.l 
[1998; Witt, Thronson & Capuano 1992) to the screening of the 
starlight by external dust (Indebetou w et al.ll2005h . It is also re- 
lated to gas and dust emission, i.e. the presence of hot stars can 
excite Biy and heat dust, which increases (H - K^) for younger 
clusters. 

The 2d index provides a way to distinguish very young clus- 
ters/star forming complexes up to 7 Myr old from the older clus- 
ters for galaxies observed in JHK^. In addition, we showed that 
the index estimates the proportion of clusters in such age ranges. 
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Table 1. Astrophysical parameters and integrated photometry of Local Group clusters. 



Cluster/Galaxy 



Q'2000 

(h:m:s) 



1^2000 

(° : ' : ") 



log [t(yr)] log[M(Mo)] E(B - V) 



d 
(kpc) 



R 

(") 



^s a-K, 



Mk j-k. 



a-j-K, J - H a-j-H H - K, cth-k. 



NGC 869/MW 


2:19:04 


57:08:06 


7.15+0.03(1 


4.20(18 


0.58 (20) 


2.3 (20) 


100 


5.06 


0.01 


-6.95 


0.10 


0.01 


-0.14 


0.01 


0.24 


0.01 




NGC 884/MW 


2:20:54 


57:08:43 


7.15+0.03(1 


4.10(18 


0.56 (20) 


2.3 (20) 


140 


5.71 


0.01 


-6.29 


0.22 


0.01 


0.10 


0.01 


0.11 


0.01 




IC 1805/MW 


2:32:43 


61:27:24 


6.30 ± 0.20 (2 


4.20(18 


0.82 (20) 


1.9(20) 


600 


3.84 


0.02 


-7.84 


0.71 


0.02 


0.49 


0.02 


0.23 


0.02 




M42/MW 


5:35:16 


-5:23:15 


6.00 ± 0.43 (3 


3.65(18 


0.05 (20) 


0.4 (20) 


900 


0.77 


0.01 


-7.26 


1.23 


0.01 


0.26 


0.01 


0.97 


0.01 




NGC2244/MW 


6:31:55 


4:56:36 


6.30 ± 0.20 (2 


3.90(18 


0.46 (20) 


1.5 (20) 


450 


4.38 


0.02 


-6.66 


0.54 


0.02 


0.27 


0.02 


0.27 


0.02 




RCW38/MW 


8:57:02 


-47:30:41 


6.00 ± 0.43 (4 


2.86(19 


3.2(21) 


1.7(21) 


400 


3.03 


0.01 


-9.23 


2.08 


0.01 


1.00 


0.01 


1.08 


0.01 




Westerlund2/MW 


10:24:02 


-57:45:29 


6.30 ±0.1 1(2 


4.00 (18 


0.79 (20) 


1.9(20) 


400 


2.85 


0.01 


-8.82 


1.62 


0.01 


0.80 


0.01 


0.82 


0.01 




Tmmpler 14/MW 


10:43:57 


-59:32:49 


6.30 + 0.43 (5 


4.00(18 


0.52 (20) 


2.7 (20) 


350 


2.99 


0.01 


-9.34 


1.03 


0.01 


0.46 


0.01 


0.57 


0.01 




NGC3576/MW 


11:11:54 


-61:18:23 


6.00 ± 0.43 (6 


2.86(19 


3.5 (6) 


2.6 (6) 


380 


3.27 


0.01 


-10.02 


1.76 


0.02 


0.86 


0.02 


0.90 


0.02 




NGC3603/MW 


11:15:07 


-61:15:36 


6.30 ± 0.43 (7 


4.10(18 


1.34 (20) 


3.6 (20) 


370 


2.38 


0.01 


-10.87 


1.51 


0.01 


0.71 


0.01 


0.79 


0.01 




Westerlund 1/MW 


16:47:04 


-45:50:38 


6.70 ± 0.09 (8 


4.50(18 


3.0 (20) 


5.2 (20) 


160 


1.76 


0.01 


-12.87 


1.61 


0.01 


0.84 


0.01 


0.77 


0.01 




[DBS2003] 179/MW 


17:11:32 


-39:10:47 


6.54 + 0.19(2 


3.80(18 


5.3 (22) 


7.9 (22) 


60 


5.47 


0.01 


-10.86 


2.10 


0.01 


1.36 


0.02 


0.74 


0.01 




Arches/MW 


17:40:50 


-28:49:19 


6.30 ± 0.20 (9 


4.30(18 


8.1 (23) 


7.6 (23) 


20 


6.40 


0.02 


-10.82 


4.35 


0.04 


2.61 


0.04 


1.74 


0.04 




Quintuplet/MW 


17:46:15 


-28:49:39 


6.60 ± 0.05 (9 


4.00(18 


8.9 (24) 


8.0 (24) 


30 


4.80 


0.01 


-12.81 


4.89 


0.02 


2.94 


0.02 


1.95 


0.02 




NGC 661 1/MW 


18:18:40 


-13:46:42 


6.48 ± 0.29 (2 


4.40 (18 


0.78 (20) 


1.8 (20) 


530 


2.05 


0.01 


-9.50 


1.94 


0.02 


1.07 


0.02 


0.87 


0.02 




M17/MW 


18:20:30 


-16:10:45 


6.30 + 0.43(10 


2.95 (19 


1.51 (20) 


1.3 (20) 


450 


1.42 


0.01 


-9.67 


1.91 


0.01 


0.83 


0.01 


1.07 


0.02 




RSGCOl/MW 


18:37:58 


-6:53:00 


7.08 + 0.06(11 


4.50(18 


7.9 (25) 


5.8 (25) 


130 


2.97 


0.01 


-13.59 


3.71 


0.01 


2.43 


0.01 


1.28 


0.01 


Q 


RSGC02/MW 


18:39:20 


-6:01:42 


7.24 + 0.09(11 


4.60(18 


4.2 (26) 


5.8 (26) 


130 


2.83 


0.01 


-12.45 


2.05 


0.01 


1.13 


0.01 


0.93 


0.01 




RSGC03/MW 


18:46:19 


-3:23:16 


7.25 + 0.05(12 


4.50(18 


4.2(12) 


6.0(12) 


200 


2.71 


0.01 


-12.64 


2.46 


0.01 


1.52 


0.01 


0.94 


0.01 




CygOB2/MW 


20:33:15 


41:18:50 


6.40 + 0.26 (2 


4.40 (18 


2.9 (27) 


1.7 (27) 


440 


2.14 


0.01 


-10.02 


1.04 


0.01 


0.60 


0.01 


0.44 


0.01 


P 


NGC7380/MW 


22:47:15 


58:06:34 


6.30 + 0.43 (2 


3.80(18 


0.60 (20) 


2.2 (20) 


300 


5.15 


0.02 


-6.77 


0.93 


0.02 


0.58 


0.02 


0.35 


0.02 




NGC330/SMC 


0:56:19 


-72:27:49 


7.40!«-2« (13 


4.56(18 


0.03 (28) 


60.0 (28) 


60 


7.24 


0.01 


-11.66 


0.74 


0.01 


0.59 


0.01 


0.15 


0.01 




NGC346/SMC 


0:59:05 


-72:10:40 


6.48 + 0.29(14 


5.60(18 


0.03 (28) 


60.0 (28) 


90 


8.23 


0.03 


-10.67 


0.63 


0.03 


0.17 


0.03 


0.46 


0.04 


NGC1711/LMC 


4:50:37 


-69:59:00 


7.70 + 0.05(15 


4.24 (18 


0.07 (28) 


50.0 (28) 


70 


8.14 


0.02 


-10.38 


0.55 


0.02 


0.43 


0.01 


0.12 


0.02 


B' 


NGC 1805/LMC 


5:02:22 


-66:06:45 


7.00!«:| (15 
7.40t|» (15 


3.52(15 


0.07 (28) 


50.0 (28) 


50 


7.52 


0.01 


-11.00 


1.01 


0.01 


0.67 


0.01 


0.34 


0.01 




NGC 1818/LMC 


5:04:14 


-66:26:02 


4.42 (18 


0.07 (28) 


50.0 (28) 


90 


7.08 


0.01 


-11.44 


0.83 


0.01 


0.63 


0.01 


0.19 


0.01 




NGC 1847/LMC 


5:07:08 


-68:58:23 


7.42 + 0.30(15 


4.44 (18 


0.07 (28) 


50.0 (28) 


40 


8.38 


0.01 


-10.14 


1.03 


0.02 


0.72 


0.02 


0.31 


0.02 


(x 


NGC 1850/LMC 


5:08:45 


-68:45:41 


7.50 + 0.20(15 


4.86(18 


0.07 (28) 


50.0 (28) 


110 


6.54 


0.01 


-11.98 


0.65 


0.01 


0.52 


0.01 


0.14 


0.01 


3 

(X 


NGC 1984/LMC 


5:27:40 


-69:08:06 


7.06 + 0.30(15 


3.38(15 


0.07 (28) 


50.0 (28) 


40 


7.90 


0.01 


-10.62 


0.74 


0.01 


0.49 


0.01 


0.24 


0.01 





Table 1. continued. 



Cluster/Galaxy 


Q'2000 

(h:m:s) 


•52000 

C:': ") 


log [t{yr)] 


log[M(Mo)] 


E{B - V) 


d 
(kpc) 


R 

(") 


ATs 


O-K, 


Mk 


J-K, 


<^J-K, 


J-H O-j^H 


H- K, 


(^H-K, 




NGC2004/LMC 


5:30:40 


-67:17:14 


7.30 + 0.20(15) 


4.36(18) 


0.07 (28) 


50.0 (28) 


100 


6.16 


0.01 


-12.36 


0.94 


0.01 


0.69 0.01 


0.25 


0.01 




NGC2011/LMC 


5:32:19 


-67:31:18 


6.99 + 0.30(15) 


3.47 (15) 


0.07 (28) 


50.0 (28) 


110 


6.60 


0.01 


-11.92 


0.96 


0.01 


0.70 0.01 


0.26 


0.01 




NGC2070/LMC 


5:38:42 


-69:06:03 


6.48!»]^ (15) 


4.78 (18) 


0.07 (28) 


50.0 (28) 


270 


4.59 


0.01 


-13.93 


0.84 


0.01 


0.31 0.01 


0.53 


0.01 




NGC2100/LMC 


5:42:08 


-69:12:38 


7.20 ±0.20 (15) 


4.36(18) 


0.07 (28) 


50.0 (28) 


170 


5.40 


0.01 


-13.12 


1.04 


0.01 


0.75 0.01 


0.29 


0.01 




NGC2136/LMC 


5:52:58 


-69:29:30 


8.00 + 0.10(15) 


4.30(18) 


0.07 (28) 


50.0 (28) 


100 


7.54 


0.02 


-10.98 


0.61 


0.02 


0.42 0.02 


0.18 


0.02 




NGC2157/LMC 


5:57:35 


-69:11:46 


7.60 + 0.20(15) 


4.31 (18) 


0.07 (28) 


50.0 (28) 


60 


8.22 


0.01 


-10.30 


0.63 


0.01 


0.52 0.01 


0.12 


0.02 




NGC2156/LMC 


5:57:50 


-68:27:38 


7.60 + 0.20(15) 


3.63(15) 


0.07 (28) 


50.0 (28) 


60 


10.31 


0.09 


-8.21 


0.24 


0.09 


0.07 0.07 


0.18 


0.11 




NGC2159/LMC 


5:58:03 


-68:37:22 


7.60 + 0.20(15) 


3.65 (15) 


0.07 (28) 


50.0 (28) 


40 


9.78 


0.03 


-8.74 


0.56 


0.03 


0.53 0.02 


0.03 


0.04 




NGC2164/LMC 


5:58:55 


-68:30:57 


7.70 + 0.20(15) 


4.18(18) 


0.07 (28) 


50.0 (28) 


90 


8.54 


0.04 


-9.98 


0.55 


0.04 


0.50 0.03 


0.04 


0.04 




NGC2172/LMC 


6:00:05 


-68:38:13 


7.60 + 0.20(15) 


3.52(15) 


0.07 (28) 


50.0 (28) 


30 


10.61 


0.04 


-7.91 


0.44 


0.05 


0.39 0.03 


0.05 


0.05 




NGC2214/LMC 


6:12:58 


-68:15:36 


7.60 + 0.20(15) 


4.03 (18) 


0.07 (28) 


50.0 (28) 


110 


7.71 


0.02 


-10.81 


0.82 


0.03 


0.57 0.02 


0.24 


0.03 




VdB0/M31 


0:40:29 


40:36:15 


7.40 + 0.30(16) 


4.85 (18) 


0.05 (28) 


783.0 (28) 


20 


12.65 


0.18 


-11.89 


0.46 


0.19 


0.46 0.14 


0.00 


0.22 




NGC604/M33 


1:34:33 


30:47:00 


6.54 + 0.06(17) 


5.00(18) 


0.04 (28) 


883.0 (28) 


40 


10.34 


0.04 


-14.40 


1.16 


0.05 


0.44 0.06 


0.72 


0.07 




Notes. Literature sources for a; 


ge, photometr 
1 (|2010|); (2) 


ic mass, reddening ; 


md distance are 


indicated within parentheses. 




|200^; 












jOi; (7) 
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